Spin echo in spinor dipolar Bose— Einstein condensates 
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We theoretically propose and numerically realize spin echo in a spinor Bose-Einstein condensate 
(BEC). We investigate the influence on the spin echo of phase separation of the condensate. The 
equation of motion of the spin density exhibits two relaxation times. We use two methods to 
separate the relaxation times and hence demonstrate a technique to reveal magnetic dipole-dipole 
interactions in spinor BECs. 
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NMR and ESR have revealed states of magnetic spin 
systems [1]. Especially the spectrum broadening contains 
informations of spin relaxation given by spin-orbit in- 
teraction, spin-spin interaction and inhomogeneous mag- 
netic fields. Usually the inhomogeneity hides the other 
characteristic effects in the system. Spin echo technique, 
however, is able to remove effects of the inhomogeneous 
diffusing. Thus the technique has been used in experi- 
ments of the magnetic spin system. 

The discovery of spin echo in 1950 by Hahn was mo- 
mentous in the field of magnetic resonance Q. Hahn 
observed the recovery of a free induction decay signal of 
spins by using the combination of three 7r/2 pulses, nam- 
ing the phenomenon "spin echo" . Carr and Purcell later 
proposed other methods of producing spin echo. The 
first, called "Method A", uses the 7r/2-7r pulse sequence, 
while in "Method B" the echoes are obtained by apply- 
ing many tt pulses after a 7r/2 pulse [3]. Spin echo has 
contributed to the development of magnetic resonance, 
especially as a method of measuring the relaxation time 
for magnetic resonance imaging (MRI). MRI is used to 
visualize the structure of human bodies for medical exam- 
inations in hospitals. Spin echo is thus a field of physics 
with practical applications. 

In low temperature physics, particularly for superfluid 
^He, NMR has made significant contributions in deter- 
mininff superfluid phases and revealing vortices and tex- 
tures Q. For example, spin super current [1,0] and Bose- 
Einstein condensation of magnons[7] have been still un- 
der investigation. Spin echo also was found in A and B 
phases in the superfluid ^He [8]. Thus magnetic reso- 
nance in superfluid ^He has been developing. 

Recently, BECs with a magnetic dipole-dipole in- 
teraction (MDDI) have received considerable attention 
0, HO, El, lill, m. The MDDI may be closely related to 
the relaxation mechanism of the spins. ^^Cr BECs have 
been realized by Griesmaier et al. p^]. A ^^Cr atom has 
a magnetic moment six times larger than that of an alkali 
atom. Hence, a ^^Cr condensate clearly shows the effects 
of the MDDI as an anisotropic free expansion which de- 
pends on the orientation of the atomic dipole moments 
[l5|. The magnitude of the MDDI can be controlled by 
modulating the s-wave scattering length [iG^. Although 



the MDDI of alkali atomic BECs is quite small, it would 
be possible to observe the MDDI in the systems as well 
as in a ^^Cr BEC []J. 

In this letter, we propose spin echo in a spinor BEC for 
the first time, studying effect of Stern-Gelrach separation 
and of MDDI by using Carr-Purcell-Meiboom-Gill se- 
quence. Investigating spin echo will provide information 
on the MDDI, which is exhibited even by ^'^Rb BECs. 

First, we briefly review relaxation in magnetic reso- 
nance and spin echo. A spin S with gyro-magnetic ratio 
7 precesses with Larmor frequency uol = jHq under a ho- 
mogeneous magnetic field H = Hqz^ which is described 
by the equation of motion of the spin dS/dt = ^[8 x H]. 
In many body systems of spin, if spins do not have an 
interaction with others, they precess with same Larmor 
frequency around a homogeneous magnetic field. On the 
other hand. The interaction between spins would make 
gradually the inhomogeneous precession through the lo- 
cal magnetic field. The decay signal, therefore, contains 
informations about spin-spin and spin-orbit interactions. 
However, if the magnetic field is inhomogeneous, the de- 
cay signal depends on also the inhomogeneity. Therefore 
it is difficult to choose effects of the interaction from the 
diffusion. Spin echo technique can remove the effect of 
inhomogeneity. 

We consider that spins without interaction are polar- 
ized to the z axis along an inhomogeneous magnetic field 
whose the spatial magnitude is almost distributed around 
Hq. Then by applying a 7r/2 pulse resonant with jHq, 
the spins tilt to the y axis with precessing. After the 
pulse, the spins precess on the x — y plane, losing the 
coherence gradually. Subsequently applying a tt pulse 
at r after the first pulse reverses the direction of spins, 
recovering the coherence gradually. The coherence recov- 
ers perfectly at r after the second pulse. This method is 
called "spin echo". On the other hand, the procedure 
can not give a perfect recovery when there is interaction 
between spins. The imperfect recovery comes from not 
the inhomogeneous magnetic field but the interaction. 

To obtain spin echo in a ^''Rb BEC trapped by a 
harmonic potential of frequency uj_\_^ we study the two- 
dimensional spin-1 Gross-Pitaevskii (CP) equations with 
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FIG. 1: (color online) Time development of the spin density vectors S projected on the x-y plane in the spin echo with 
(jj±r = 5. The magnitude of S increases from blue to red. 



the MDDI: 

+-^^ / dr'^-f^S,{r')Si,^,. (1) 

Here Vtrap = Mc<;^(x^ -\- y'^)/2 is the trapping potential, 
/i the chemical potential, g the Lande's g-factor, jub the 
Bohr magneton, and 

5i=EV':^a/3# (2) 

a,f3 

is the component of the spin density vector S repre- 
sented by the spin matrices 5^^. The short range in- 
teraction constants are cq = 47rfi^(ao + 2a2)/3M and 
C2 = 47rfi^(a2 — ao)/3M given by the s-wave scattering 
lengths ao and a2. The long range interaction constant 
is Cdd = Mo^^Mb/^'^ with vacuum magnetic permeabil- 
ity /io- For ^''Rb, the parameters satisfy the relation 
Cdd < — C2 <C Co [12|. We investigate the spin echo under 
a gradient magnetic field H = {Gx -\- Ho)z with constant 
G. 

The dynamics start from the stationary states with 
\cdd/c2\ = and 0.02 in a gradient magnetic field 
of magnitude 0.98i^o < |H| < 1.02i^o. Applying a 
Urot = Hi{cos{jHQt)x — sm{'y HQt)y} with Hi = O.IHq, 
we calculate the GP equation using the Crank-Nicolson 
method, obtaining the spin echoes and effects of dephas- 
ing in the echo. 

Spin echo is clearly observed in Fig. 1, showing the 
dynamics of S projected on a x-y plane for Cdd = and 
uj±r = 5. After the 7r/2 pulse, the spins are oriented 
in the -y axis (Fig. [IJa)) and start precessing on the 
plane. Then, the precessions gradually diffuse because 
the Larmor frequency is spatially dependent due to the 
gradient magnetic field (Fig. HKb), (c)). At t = r + 



t^/2 (Fig. md)), the tt pulse is applied, which reverses 
the direction of the spin (Fig. [D^e)). The spins then 
gradually become coherent (Fig. [D^f), (g)), eventually 
refocusing at tpeak = 2r + ^7^/2 +^7r(Fig. [Uh)). After the 
echo, the spins start to defocus again. These features are 
represented in the time development of the expectation 
value {Sy^,) = J dri/j^^Sl'^i/jp in Fig. [2 The signal of 

(Sy) gradually decays from t^/2 (a) to r + t^/2 (d) and 
then increases from r + t^/2 + tn (e) to 2r + 3t^/2 (h). 

To show that this is genuine spin echo, we investigate 
tpeak by changing r, thus confirming the spin echo for 
\cdd/c2\ = and 0.02. Figure [3] indicates that the re- 
lation is satisfied until ujj_r = 6.5 for Cdd = and 4.5 
for \cdd/c2\ = 0.02. These deviations from the relation 
tpeak = 2r + 3^7^/2 cau be understood by investigating the 
equation of motion of spin based on the GP equations. 
From Eqs. ([T]) and ([2|), we obtain the equation of motion 

^ = V • J. + 7 [S X H], + '-flV,,S,^^,,^, (3) 

where n/(2Mz)S'^^(V^^VV^« - is term of 

spin current, which is similar to spin current in super- 

fulid ^He [18], and = J dr' ^^^^^^Sk{r'). The first, 
second, and third terms of Eq. (|3]) are derived from the 
kinetic, Zeeman, and MDDI terms of the GP equation. 
The differences between the circles and triangles in Fig. 
[3] are due to the dipole term, while their deviations from 
tpeak = 2r + 3^7^/2 ^rc caused by the spin current given by 
separation of the condensates due to the gradient mag- 
netic field. 

The magnetic resonance of this system is greatly af- 
fected by separation of spinor condensates. As shown in 
Fig. m (a), (b), (c), the condensates separate under the 
gradient field after the 7r/2 pulse, keeping approximate 
symmetry between x > and x < 0. Applying the tt 
pulse changes the ipi and i/j-i components to ip-i and 
ipi respectively, which reverses S. The tt pulse breaks the 
symmetry of the distribution of the condensates. After 
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FIG. 2: Time development of {Sy)/h (black line) and {Sz)/h 
(white circles). The symbols (a)-(h) indicate the times shown 
in Fig. [1] The gray zones represent the intervals of the pulses. 
The behavior of {Sz)/h shows that the Ti relaxation is irrel- 
evant to the dynamics {Sy)/h. 



the pulse, ijj±i move to the center, tending to overlap 
each other for u±r < 7.5. The overlap is not clear for 
uj±T > 7.5. The effect is caused by a repulsive inter- 
action between the different components, which is eas- 
ily found in the GP equations of Eq. ([T|). Thus, spin 
echo works properly for pulses shorter than u;±r = 7.5, 
because phase separation is not yet significant. On the 
other hand, the MDDI operates as an attraction between 
V^i and ip-i and makes the echo peaks appear earlier than 
the echo without the MDDI, as shown in Fig. [3l 

We now apply Methods A and B, which are major se- 
quences in NMR and ESR for measuring two relaxation 
times, to spinor condensates. The dynamics of Eq. ([3|) 
is characterized by two relaxation times. One is the re- 
laxation time Tnu due to the non-uniform condensates 
and the other is the dipolar relaxation time T^^. In 
Method A, the decay of echo peaks is measured using 
a 7r/2-7r pulse sequence for various values of r. In liquids 
and gasses, the decay is generally caused by atomic self- 
diffusion; the first term of Eq. (|3]) operates similarly in 
our systems. In Method B, a 7r/2-7r-7r . . . pulse sequence 
is used, with n times tt pulses applied after a 7r/2 pulse. 
Carr et al introduced the effect of diffusion on liquids 
and gasses into the decay of a signal by using a random 
walk model and obtained transversal decay of the 7r/2-7r 
sequence exp[— t/T2 — 7^a^G^t^/24n^r] [3], where a is a 
fixed distance of the random walk. Thus, the diffusion 
decreases with increasing n. The diffusion is given by 
the first term of Eq. ([3]) in our case. For investigat- 
ing Tnu and Tdd^ we performed a numerical simulation 
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FIG. 3: Time oo±tpeak at which echoes appear versus the 
interval co±t between 7v/2 and tt pulses. The circles and tri- 
angles are the results for \cdd/c2 \ = and 0.02, respectively. 



of Methods A and B, shown in Fig. [5l Figure [5] (a) 
shows multiple exposures of {Sy)/h at t = tpeak for vari- 
ous T for Method A. The echo peaks decay with r. When 
the condensates separate for large r, the decay does not 
satisfy the exponential behavior. Thus, we fit the ex- 
ponential formula to the echo peaks until uj±r = 7.5, 
obtaining a diffusion time of uj±t = 40 when {Sy)/h de- 
cays to e~^. This diffusion time should be equal to Tnu- 
Figure [5] (b) shows a single exposure of the echo peaks for 
Method B with uj±r = 3.5 and n = 17. In the simulation, 
rather than applying the original 7r/2-7r-7r. . . sequence, 
we apply a 7r/2-7r-(— 7r)-7r-(— tt) . . . sequence, considered 
by Meiboom and Gill [l9[, to cancel deviation due to re- 
peated TT pulses. In the case of no MDDI, the echo peaks 
do not decay, showing that the Method B sequence re- 
moves the effect of phase separation. While the signal 
with the MDDI decays, the relaxation time should be 
given by uoxTdd ~ 500. These results show Tnu < Tdd- 
Method B enables us to reveal the MDDI through Tdd in 
spinor dipolar BECs even if the MDDI is too weak to be 
easily detected. 

In conclusion, we have numerically realized spin echo 
in a trapped ^''Rb BEG by calculating the spin-1 two- 
dimensional GP equation. We have investigated how the 
spin echoes are affected by the phase separation and the 
dipole-dipole interaction. The equation of motion of the 
spin density (Eq. (|3|)) derived from the GP equation 
shows two relaxation times Tnu and Tdd- In order to dis- 
tinguish Tdd from Tnu by using spin echo, we introduce 
Methods A and B, obtaining the relation Tnu > Tdd- 
The study of spin echo of spinor BEGs will help develop 
the field of magnetic resonance in cold atomic BEGs sys- 
tem and may be applicable to MRI in these systems. 



4 



(jj±T = 5 (jj±r = 7.5 cjx^ = 12.5 




X X 3; 



FIG. 4: Density profiles \ijjc,{y = 0)|^ = aI|V^c.(2/ = 0)|V^ 
witii ah = ^JW[mZjx_ and x = x/a^ without the MDDI just 
before the tt pulse ((a), (b), (c)) and just after the tt pulse 
((d), (e), (f)) for cj^r = 5, 7.5, and 12.5. 




ijj\t (jj±t 
(a) Method A (b) Method B 

FIG. 5: (a) Results of Method A, showing multiple 
exposures of {Sy)/h with (circles) and without (trian- 
gles) the MDDI at tpeak for various r, where ix)_\_T — 
3.5, 4.0, . . . , 7.5, 10, 12.5, 15, 17.5. The black line shows the 
Carr-Purceh form exp[-6t - c^t^] with h/c = 0.072. (b) Re- 
sults of Method B, showing a single exposure of the signal of 
peaks with cj^r = 3.5 and n — 17. The line exp(— 6t) is fit to 
the early peaks, which are little affected by exp(— c'^t^/17^). 
We obtain ijOi_Tdd ^ 500. 



Although this work is limited to the two-dimensional 
case, three-dimensional spinor BECs, especially with the 
MDDI, should show interesting features depending on ge- 



ometry. A study of three-dimensional spinor BECs will 
be reported soon elsewhere. We also would like to study 
spin wave and spin transport as analogy to magnetic res- 
onance of superfluid ^He. 
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